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Abstract 

In order to derive a large set of Hamiltonian dynamical systems, 
but with only first order Lagrangian, we resort to the formulation in 
terms of Lagrange-Souriau 2-form formalism. A wide class of sys- 
tems derived in different phenomenological contexts are covered. The 
non-commutativity of the particle position coordinates are a natural 
consequence. Some explicit examples are considered. 



* e-mail: Luigi.Martina@le.infn.it 



1 



1 Introduction 



Recently interest in dynamical systems with non commuting (not necessarily 
non canonical) variables stems from Condensed Matter Physics [JJ, Optics 
[2] and String Theory [3]. The applications include a Bloch electrons [3], 
the Anomalous Hall Effect [5] , the Spin Hall Effect [6] , and the Optical Hall 
effect [?]. For example, 

a) The semiclassical equations of motion in the e n (k) energy band of a 
Bloch electron in a crystal solid read 

P= ^"ffl -k x 6(fc), k = -eE -efx B{r). (1.1) 
dk 

The purely momentum-dependent @i(k) = e^id? Ai(k) is the curva- 
ture associated to the so-called Berry connection A. 

b) The semiclassical equations that describe the spin-Hall effect into semi- 
conductors near the degenerate point k = of the valence band read 

f=^^- + %x@ s , %=-eE, (1.2) 
dk 

where s = ±| is the spin helicity of the holes, the energy is E s (k) = 
(A — Bs 2 ) k 2 and the Berry curvature due to the lattice structure 

is & s = s (2s 2 — |) A. The trajectories followed by opposite helicity 
holes separate during the motion, providing a tool for spintronic. 

c) The optical Magnus and the optical Hall effects [21 [71 E] are described 
by the approximate equations 

si- 1 

r m p grad(— ) x p, p w — n 3 w 2 grad(— ), (1-3) 

co n n 

where s denotes the photon's spin, parametrizing a term which devi- 
ates the light's trajectory from the predictions of ordinary geometrical 
optics. 

d) The motion of a Bogoliubov quasiparticle in a superfluid vortex [9] is 
described by 

dh dh 
Mq + F xf= - — , Mr=— , (1.4) 
or oq 
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where the effective mass Mij = 5ij — and magnetic field F{j = 

f^f — are defined in terms of a gauge potential A = A(q,r), and 
h = H (q,r) is some Hamiltonian. 

e) The string theory inspired dynamical system [10] 

„ = « + e„^, B= *I + me ,^, (i.5) 

m oxj oxj oxiOXj 

for the Kepler potential V oc r _1 in a weak non-commutative back- 
ground predicts a perihelion point precession in the planetary mo- 
tions, providing us with a measurable quantity to test cosmological 
models. 

In all these models, the momentum satisfies a first order ODE in terms of a 
(possibly of the Lorentz-type) force. But the same structure appears in the 
equations for the position, which are distinguished by an unusual velocity 
relation, as a result from the Berry phase contributions, associated with the 
environment in which the particle moves, or by a postulated fundamental 
area scale. Thus, the momentum and the velocity may not be proportional. 

In the absence of any magnetic field, for example, equations ([I]) reduce 
to those used to explain the AHE observed in ferric materials [5j: the effect 
comes, entirely, from the anomalous velocity term. Let us emphasize that 
such a behavior has been advocated a long time ago [11] . It has been shown 
indeed that the Lorentz-equation follows from general principles, but the 
momentum- velocity relation is model-dependent. Relations like p = mr are 
indeed mere Ansatze, and it is not required by any first principle. 

More recently, similar features has been found in 2D systems for the so 
called exotic particles, i.e., ones which are associated with the kinematical 
two-parameter central extension of the planar Galilei group [121 113] . These 
models are, once again, instrumental in describing the Fractional Quantum 
Hall Effect. Those systems support a sort of duality between the magnetic 
field, B(f), and its analog in momentum space, @{p), called sometimes a 
dual magnetic field. By consequence the Poisson bracket of the position 
variables they may no longer vanish, in analogy with the components of the 
momentum in the presence of a magnetic field. In the relativistic case, sim- 
ilar, but more elaborated constructions have been proposed [144 115] . which 
however we do not discuss here. 

Even the question of consistency of the above effective models with the 
general principles of mechanics is legitimate, since they are mainly derived 
by some semiclassical de-quantization procedure, which does not necessarily 
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fit into the framework of classical mechanics. Thus, our primary goal below 
is to prove that no new mechanics has to be invented : all these models 
fit perfectly into Souriau's presymplectic framework of Classical Mechanics 
|16j . The second aim is to formulate the hamiltonian theory in that con- 
text, exploiting the supplementary structure encoded into the second central 
extension of the Galilei group. 



2 The Lagrange-Souriau 2-form 

The modern geometrical formulation [16] of the calculus of variations, origi- 
nated by Lagrange and continued by Cartan [17], consists in mapping the La- 
grangian function L = L(x, v, t), defined on the evolution space TQxM — > M, 
into the so-called Cartan 1-form A, defined by the crucial property 

L( 7 (t), 7, t)dt = J A with A = — dx i + - — Pi j dt. (2.1) 

where 7 = (7(4), 7, i) is the lifted world-line in the tangent bundle of the 
evolution space TQ x R. The exterior derivative of the Cartan form provides 
us with a closed Lagrange-Souriau (LS) 2-form a = dX, which in general 
cannot be separated canonically into a symplectic and a Hamiltonian part 
[16] . However, the associated Euler-Lagrange equations are still expressed 
by determining the kernel of a, i.e. by the equation (7(7, •) = 0. If the kernel 
is one-dimensional, the variational problem is called regular. Otherwise 
the singular case requires to resort to the symplectic reduction techniques 
[16l \T7\ \18\ [T9] . Thus more general procedures have to be adopted to build 
such a system and clarify their Hamiltonian structure. 

Conversely, again following Souriau [16] , a generalized mechanical system 
is given by a 2-form a, which is closed da = and with constant rank d. 
Then, its kernel defines an integrable foliation with (i-dimensional leaves, 
which can be viewed as generalized solutions of the variational problem. 
Moreover, by the Poincare lemma, da = implies the local existence of a 
Cartan 1-form A. Rewriting it as A = a a d^ a , one can plainly define the 
first-order Lagrangian function on the evolution space as 

£ = a a £ a such that j Cdt = j A. (2.2) 

Thus, the above-mentioned models do not have a usual Lagrange function 
defined on the tangent bundle. Put in another way, the position does not 
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satisfy a second-order Newton equation. The general question of the exis- 
tence of a locally defined Lagrangian has been discussed [20j EH] ■ 
Moreover, if the LS 2-form a can be split as 

a = u - dH A dt, (2.3) 

where w is a closed and regular 2-form on the phase space TQ and H is a 
Hamiltonian function on TQ x R, than the equations of motion read 

u (7) = dH. (2.4) 

Assuming that u is regular, the inverse (ui a P) of the symplectic matrix oj = 
io a /3 ( i.e. ui al3 LUpj = 5") yields the Hamilton equations £j = H}, through 
the Poisson brackets defined by 

{f,g} = u aft d a fd p g. (2.5) 

In the case of singularity, as said above, symplectic reductions are needed. 

Thus, in Souriau's framework, one can state the inverse problem of the 
calculus of variations for given equations for a one-particle system in pres- 
ence of a position-dependent force field E only, rewriting them as the set of 
1-forms on TQ x R 

ai = df — p dt ai= vn dp — E dt, (2-6) 

and view both of them as hyperplanes in the evolution space, defined by 
the kernel of the one- forms ctj. Notice that we have introduced the kinetic 
momentum p = mv, and we will refer to it in the sequel. Then, the si- 
multaneous solutions correspond to the intersection of these hyperplanes, 
described by the kernel a(5y,-) = of exterior product a = a\ A 02, 
where A warrants the antisymmetry. In presence of an electromagnetic field 
E(f,p,t) ,B(f,p,t) acting on a charge e, Souriau [16] generalized the pre- 
vious simplest 2-form to 

a = (mdv - eEdt^j A (df - vdt) + eB ■ df x df, (2.7) 

where we have defined (dfx df) k = ^kijdxi A dxj. Then, the usual equa- 
tions of motion of a charged particle in the electromagnetic field are seen 
to arise as the kernel of a, together with the regularity condition da = 
leading to the the homogeneous Maxwell equations. These formulas can be 
readily generalized to the multi-particles case. 

Now, in the same spirit we would like to write down a Lagrangian 2-form 
for a particle of mass m, which is subjected both to the electromagnetic field 
E (f, t) and B (f, t), and to a peculiar environment, the local characteristics 
of which may depend on the momentum. 
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3 The (2+l)-dim Duval - Horvathy model 



The simplest example of such a situation is the exotic mechanical model in 
2 space dimensions proposed in [13], defined by a form generalizing (|2.7p 
which can be split as in ()2.3p : precisely we introduce 

1 f P 2 \ 
a = dpi A dxi + -9 €ij dpi A dpj + eB dx\ A dx2 + d I V eV I A dt, (3.1) 

where B (r) is the magnetic field, perpendicular to the plane, V (r) the elec- 
tric potential (both of them assumed to be time-independent for simplicity) 
and 9 is a constant, called the non- commutative parameter, for reasons clar- 
ified below. The resulting equations of motion read |13j 



m*Xi = pi — em9 tijEj, p\ = eEi + eB €ijij, (3-2) 

where we have introduced the effective mass m* = m (1 — e 9 B). The novel 
physical features are: i) the anomalous velocity term —em9eijEj, so that 
Xi and pi are not in general parallel, ii) the interplay between the exotic 
structure and the magnetic field, yielding the effective mass m* in ([3]). 

The LS 2-form (|3.ip can obtained by exterior derivation from a Cartan 
1-form on the evolution space M 2 x M, namely from 

p 2 9 

A = (Pi- Ai )dxi - — dt + - €ij pi dpj (3.3) 

defining the action functional as in (|2.2 j) . Accordingly, for m* ^ 0, the 
associated Poisson brackets are 

m Til m 

{x 1 ,x 2 } = — 9, {xi,pj} = —5ij, {pi,p 2 } = —eB. (3.4) 
m* m* m* 

and they automatically satisfy the Jacobi identity. When effective mass van- 
ishes, i.e. when the magnetic field takes the critical value B cr n = \, the 
system becomes singular. Then symplectic reduction procedure leads to a 
two-dimensional system characterized by the remarkable Poisson structure 
{21,22} = 9, reminiscent of the "Chern-Simons mechanics" |22| . Thus, the 
symplectic plane plays, simultaneously, the role of both configuration and 
phase spaces. The only motions are those following the Hall law. Moreover, 
in the quantization of the reduced system, not only the position operators no 
longer commute, but the quantized equation of motions yields the Laughlin 
wave functions [23], which are the ground states in the Fractional Quan- 
tum Hall Effect (FQHE). Thus one can claim that the classical counterpart 
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of the anyons are in fact the exotic particles in the system (|3,2p . In the 
review article [23] several examples of 2-dimensional models, which gener- 
alize the form f)3. 1 j) and the equations (|3.2p have been discussed. Here let 
us conclude that the Poisson structure (|3.4p can be obtained by applying 
the Lie-algebraic Kirillov-Kostant-Souriau method for constructing dynam- 
ical systems to the (2+1) 2-fold centrally extended Galilei group |25t 126] . 
The two cohomological parameters are the mass and a second, "exotic", 
parameter, identified with a constant Berry curvature in the context of the 
condensed matter physics, or as a non relativistic limit of spin (see [23] and 
references therein). The latter is highlighted by the non-commutativity of 
Galilean boost generators, 

[K 1 ,K 2 ] =i6m 2 . (3.5) 



4 The general model in (3+l)-dim 

After the above digression on two-dimensional models, let us look for further 
generalizations [2Z] - [31] of the 2-form (|2.7p with momentum dependent 
fields. Straightforward algebraic considerations lead to define on the space - 
momentum - time evolution space the manifestly anti-symmetric covariant 
2-tensor on the evolution space of 

a = [(1 - m) dpi - e Ei alt] A (dr« - gi dt) + f e B k e kij dr { A drj + 

3 Kfc ekij dpi A dpj + q k e kij dr { A dpj, (4.1) 

where we have put into evidence the Lorentz contribution to the Lagrangian, 
the quantities g, k and q are 3- vectors and Q = diag (/Xj) is a 3 x 3 matrix, 
respectively. They depend on all independent variables and have to be de- 
termined in such a way that a is closed and has constant rank. By the 
expression 1 — m we would like to distinguish among the bare mass, normal- 
ized to 1, and the possible variable contributions. It is interesting to note 
that, with respect to the expressions adopted in [16] and [13], we have to 
introduce the 1-form dri — gi dt, which defines a new conjugate momentum. 

The equations of motion can be written as the kernel of the LS 2-form 
a {Sy, •) = for a tangent vector Sy = (5r,5p,5t). Specifically one obtains 
the equations 

e8r-E-(l-Q)8p-g = 0, (4.2) 
(l-Q-Q A )8p = e(E5t+5rxB), (4.3) 
(1 - Q) (Sf-g St) = -5fxq-5pxK, (4.4) 
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where we have introduced the anti-symmetric matrix (QA)ij = e ijk Qk- From 
the equation (|4.3|) we can solve with the respect of the vector 5p, if the mass 
matrix M = 1 — Q — Qa is invertible, i.e. if det (M) 7^ 0. Under such an 
hypothesis together with det (m — Q) 7^ 0, one replaces dp in the equation 
), finding an equation for the position tangent vector Sr of the form 



M* Sr= (Jl-Q)g + e KM ■ Ej St, (4.5) 

where the effective mass matrix is given by 

M* = M+ (2Q A - e QM~ l B) , Q {j = e ijk K k , B {j = e ijk B k . (4.6) 

Both M* and the eq. (|4,5|) generalize of the expressions obtained in [13] 
leading to (13. 2|) . Singularities in the motion can arise, both from the van- 
ishing of det (M) and degeneracies of the first factor in (|4.6p . that is at 
det (M*) = 0. However, if this is not the case, one can solve (|4.5p w.r.t Sr 
and show that 



1. equation (|4.2p is identically satisfied independently from the specific 
choice for the vector g, 

2. the equation (j4.3j) becomes 



M " S ? = D^Mj{ RB -t TNB ) St ' < 4J > 

where matrices R and N have an involved dependency on m, B, k, Q 
and Qa to be spelled here. 



Thus we have a system of simultaneous first order differential equations: 
(|4.5p for the velocity of the particle 4? and (|4.7p for the momentum variation 

Notice that these two equations simplify to the equations (8) and (9) of 
the [13] , when Q and Qa go to and g = p. 

Now, the closure condition in the evolution space da = ( called also 
"Maxwell Principle" by Souriau |16| ) of the Lagrangian-Souriau form a in 
(|4.ip leads to the coefficients of the independent 3-forms. It is quite natu- 
ral to assume the following limitations on the involved functions: d pi Ej = 
d Pi Bj = . Thus we are lead to the equations 

d r .Bj = 0, e kij d n Ej = -d t B k , (4.8) 

d p .Kj = 0, e kij d n [(1 - Hj) gj] = d t K k , (4.9) 

dtm = d n [(1 - ^ gi] , \e kij d n [(1 - fij) gj] = d t q k , (4.10) 

d ri Hj = £ij k d rj q k , d r% Kj = Eij k d Pk Hi + d Pi qj — Sijd Pk q k , (4.11) 
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with the residual closure relations (without summation over repeated in- 
dices) 



a rj [(l-ii i )9i] + d ri [(l-n j )g j ] = 0, i^j = 1,2,3. (4.12) 



One can observe that the homogeneous Maxwell equations (|4.8p are the only 



analogs of the previous relations in the velocity (momentum) space for the 
vector-field k, which measures the extent to which the spatial coordinates 
fail to commute in three dimensions, as we will see in the Hamiltonian 
formalism. For such a reason, sometimes Ft is called dual magnetic field. 

If Ft is non trivial in time, then a change in the velocity dependence is 
induced for the mass flow (1 — Q)g, as prescribed by the second equation 
in (14. 9p . In its turn, equations (14.10j) say how the particle mass may change 
in time. This seems to be a quite unusual situation, but we cannot discard 
it at the moment. On the other hand, the first set of three equations in 
(|4.10p has the form of independent continuity equations, leading to the global 
conservation law for the total mass, i.e. dt (X^Mi) + ®n [(1 — f^i) 9i] = 0> 
which however holds separately in different directions. Due to q, also the 
skew-symmetric contributions to the mass matrix M may change on time, 
but they generate modification of the mass flux in space, accordingly to the 
second set of equations in (|4.10p . 

The equations in (14. lip are more difficult to interpret: they provide con- 
sistency relations for both the space and the velocity dependency among the 
mass matrix elements and the dual magnetic field. Putting such expressions 
into the equation of motion in the form (|4.3p - (|4.4p . in a pure axiomatic way 
one re-obtains the equations found in the context of the semiclassical motion 
of electronic wave-packets in [2] . 

For a particle in a flat evolution space, i.e. for ^ = and momen- 
tum pi = g, L , one easily concludes that k and q have to be constants. 
The analysis of the closure relations (|4.10|) - (|4.11[) leads to the expression 
K i = Ylj^i (xjdpjQi - XidpjQj) + Xi ,where the g^'s and \i 's depend only on 
the velocities and moreover the divergenceless condition d Pj Xj = has to 
be satisfied. A remarkable example for its phenomenological implications 
is provided by the monopole field in momentum space k = Ojj^s, which is 
indeed the only possibility consistent with the spherical symmetry and the 
canonical relations {Vj,pj} = 5{j [29J. Its expression appears to be consis- 
tent, at least qualitatively, with the data reported in ([5]) and in Spin Hall 
Effects [6], [30]. 

Limiting ourselves to two spacial dimensions and setting k% = const, a 



restrictions 




Equations (|4.9p are the 
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free particle admits in fact the "exotic" two-parameter centrally extended 
Galilei group as symmetry |13j . In a more general situation one deals with 
a momentum-dependent non-commutativity field K3 = n(p), like it was con- 
sidered by [32] with k = —j^p7^ijPi r j (and /Xj = ), or K a ^ = § ^2^3/2 f° r 
the planar relativistic model for a spinning particle (again a sort of monopole 
in relativistic momentum space) [331 02] • 

Finally, in the singular submanifold of the phase space defined by M* = 
0, we need to look at the proper restrictions on vector-fields 5p and 5r, in 
order to avoid motion with infinite velocities. In particular in two dimen- 
sions, assuming that the only non vanishing components of the magnetic, of 
the dual magnetic fields and of the q are the ones orthogonal to the plane of 
motion, with det (1 — Q) 7^ 0, the singularity mass manifold is described by 

B cnt = g2 + (m -^ l)(m -^ 2) . (4.13) 

For sake of simplicity, here we assume all the above quantities as constants 
on space- velocity variables. The equations of motion (|4.2p - (|4.4p and the 
closure relations are satisfied by the following constraints for the mass flow 
and electric field components 

(I- a) a = ^^Ej-ejjEjjl- fn)) _ (qEj - gjEj (1 - /i t )) 

I g 2 + (! _ Ml ) (x _ M2 ) Bcrit ' ^ • > 

(1 - H2) d ri E x + (m - 1) 8 r2 E 2 + q {d ri E 2 + <9 r2 £i) = 0. (4.15) 

Notice that the last relation comes from the closure condition on the mass 
terms (|4.12p . not from the Faraday law (|4.8p . which is identically satisfied. 
Since the same equations provides the time evolution of the mass terms and 
of k, the only left by above assumptions, the electric field can depends at 
most linearly on space coordinates, compatibly with (|4,15p . Finally, the 
equation (|4.14|) generalizes the Hall law discussed in the previous Section. 



5 Hamiltonian Structure 

Comparing the system (|4.5p - (|4.7p with the previous ones in (|3.2p . one rec- 
ognizes the general common structure, due to derivation from the same 
unifying Lagrange-Souriau approach provided by the 2-form (|4.ip . even if 
the peculiarity of the doubled folded central extensions is enjoyed only in 
the 2-dimensional setting. 
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The 2-form a can be obtained as the exterior derivative of the Cartan 
1-form, 

A= (p + ~£\ -dr + K -dp-Tdt. (5.1) 

In this formula the field A (r, t) is the usual electromagnetic potential, 
such that B = V f -? x A, for which we have postulated only a space - 
time dependency, eventually resorting to a suitable gauge transformation. 
Then, the electric field (in fact any space-time dependent force), is given by 
E = -V r -T - d t~2, where we assume the following decomposition for the 
scalar function T(r,p, t) = 6 (p, t) + ip (r, t), in order to keep valid the previ- 
ous restrictions on the dependency of E fields. On the other hand, the field 
1Z (r,p, t) defines the dual magnetic field k = Vp- x 7Z, the mass flow is given 
by (1 - Q)g = —VpT— d t TZ, m = d r JZi and q k = d ri TZj = -d rj TZi, where 
k, i,j are in cyclic order. The last equality imposes a set of constraints on 1Z 
in such a way the Lagrangian-Souriau 2-form takes exactly the expression 
(|4.ip . The above restrictions implies certain second order relations, which 
assure that also the first set of closure relations in (|4.1ip are satisfied, namely 

d n d r] TZ k = (i,j,k cyclic), d%.Kj = -d n d r .Hi (i ^ j) , (5.2) 

with no summation over repeated indices in the last equation. The remain- 
ing closure relations in (|4,lip are identically satisfied, while those in (|4.10p 
require the supplementary constraints on time derivatives 

d t (d n TZ 3 + 2d Pj Ki) =0 i^j. (5.3) 

Due to the special form we assumed on the force and magnetic fields, the 
above restrictions on 1Z limit its space-time dependency, leaving however the 
gauge freedom with respect the momentum variables. 

Thus, in terms of the potentials, or connections, introduced in (|5.ip . the 
equations of motion (I4.3p - (|4.4[) become 

(Sij + Eij) fj + Qij pj = d Vi £ + d t TZi, 

Bij fj + (5ij + Hy) pj = -d n (p-d t Ai, (5.4) 

where we have used the matrices (see also eq. (14. 6p ) 



I j^. i> j ^ _ £ ijk^jhkdr h Ak, &ij — £ijk£jhkdp h 7Zk(5.5) 



Differently from system (|4.3p - ([4.4p on the evolution space, now the dynam- 
ical system (]5.4p is defined on the tangent manifold of the phase space 
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TQ = {£ = (f,p)} space. But, if dtA = d{R- = 0, it is possible to rear- 
range the form (|4,ip as a = uj — dH A dt by introducing the symplectic 
2-form on TQ 

u = + Sij ) dr i A dpj + - [By dr^Adrj- 9y dp i A dp j] (5.6) 

and the Hamiltonian function H = T = £ (p,t) + (p (f,t). 

Thus, the closure of w = w aj 8 d£ Q A is assured by that one of a, i.e. 
by the equations (|4.8p - (|4.12p plus the restriction (|5.2p . Then, the space 
TQ becomes a Poisson manifold, with Poisson brackets defined, as in (|2.5|) . 
by the co-symplectic matrix 



" ; - ^l-|Tr(H 2 + X(l + 2H)G) ) 1 (5.7) 



+ [i-}n(s» + Be)](_° 1 J) + 



+ [ S ' ] ^ , n i^^2 , bq m f o l 
-S+[3,B] 

( o (H 2 + ie) T 

l-(S 2 + B9) 



non degenerate for y/deb (v a p) = 1 - ^Tr (H 2 + B (1 + 2 S) 9) ^0. Such 
a factor generalizes the denominators present in the Poisson brackets (|3.4p 
or (14. 6|) . Moreover, it crucially appears in the expression of the invariant 
phase - space volume, ensuring the validity of the Liouville theorem. Fi- 
nally, notice that the Poisson structure is determined only by gauge invari- 
ant quantities and brackets involving position coordinates rj are in general 
non-commutative. As it has been elsewhere remarked [24], it is possible 
to perform a change of variables leading to commutative position variables 
by a point transformation of the form T{ — > r • = rj — TZi(f,p). However, 
the vector field 1Z is defined up to a gauge transformation generated by an 
arbitrary function on (f,p). Thus, the meaning of the notion of position is 
unclear in such a context. 

For a charge subjected only to a monopole of strength 9 in momentum 
space and to a uniform electric field, the equations of motion obtained (|5.8p 



are readily integrated [30] . The particle suffers a shift A = |^ in the direction 

E x po, being pq the initial linear kinetic momentum. This quantify the 
discussion in point b) in the Introduction. On the other hand, if the charge of 

Hamiltonian H = Mj- is driven only by a magnetic and by a dual monopole, 
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the equations are 



M* 



Pi - e6 




(5.8) 



|f| 3 |p| 3 



\p\\f\ 



M* 



Pj r k 



(5.9) 



r| 3 |p| 3 




where M* = |r| 3 |p| 3 — e9 r ■ p. 

A final remark in connection with the coupling to the electromagnetic 
field adopted in (|4.ip . sometime said the minimal addition and leading to the 
Hamiltonian T above. It is quite different from the usual minimal coupling 
procedure and it yields a very different Poisson structure. In the context of 
the 2-dimensional systems this problem was reviewed in [23]. In particular 
the two formulations, in that context, were proved to be equivalent under a 
classical Seiberg-Witten transformation of electromagnetic fields. No results 
are yet available in the situations discussed in the present paper. 

In conclusion a wide set of dynamical systems is derived from the Lagrange- 
Souriau approach in 3-dimensions. Generalizations to higher number of de- 
grees of freedom seems straightforward. We have shown the conditions to 
assure their Hamiltonian formulation. From which an analysis for their in- 
tegrability properties can be pursued more plainly, by resorting to standard 
methods. Alternatively, one can perform a symmetry analysis directly on 
the LS 2 -form (|4"TT]) accordingly to [25] and [33]. 

Acknowledgments 

The author expresses his indebtedness to P. Horvathy for having been intro- 
duced to the problem and for continuous and stimulating discussions. Some 
of the results were obtained in collaboration with C. Duval, Z. Horvath and 
P. Stichel. The work has been partially supported by the INFN - Sezione of 
Lecce under the project LE41. 



References 

[1] D. Xiao, M. C. Chang and Q. Niu, Rev. Mod. Phys. 82 (2010) , 1959 
[2] K. Y. Bliokh, J. Opt. A: Pure Appl. Opt. 11 (2009) 094009. 
[3] R. J. Szabo, Phys. Rept. 378 (2003) 207. 



13 



[4] A. Bohm, A. Mostafazadeh, H. Koizumi, Q. Niu and J. Zwanziger, 
Springer Verlag (2003). 

[5] Z. Fang , N. Nagaosa , K. S. Takahashi ,A. Asamitsu , R. Mathieu , 
T. Ogasawara , H. Yamada , M. Kawasaki , Y. Tokura , K. Terakura , 
Science 302 (2003), 92 

[6] S. Murakami, N. Nagaosa, and S.-C. Zhang, Science 301, 1348 (2003). 

[7] K. Yu. Bliokh, Yu. P. Bliokh, JETP Lett. 79, 519 (2004); K. Yu. Bliokh 
and Yu. P. Bliokh, Phys. Lett. A333, 181 (2004) [physics/0402110]; 
N.R. Sadykov, Theor Math. Phys. 135, 685 (2003) ; M. Onoda, S. Mu- 
rakami, N. Nagaosa, Phys. Rev. Lett. 93, 083901 (2004); 

[8] C. Duval, Z. Horvath, P. A. Horvathy, Phys. Rev. D74, 021701 (2006) C. 
Duval, Z. Horvath, P. A. Horvathy, Journ. Geom. Phys. 57, 925 (2007); 

[9] C.-W. Zhang, A. M. Dudarev, Q. Niu, Phys. Rev. Lett. 97, 040401 
(2006). 

[10] J. M. Romero and J. D. Vergara, Mod. Phys. Lett. A 18 (2003) 1673 

[11] W. G. Dixon, II Nuovo Cimento 38, 1616 (1965); J.-M. Souriau, Ann. 
Inst. Henri Poincare, 20 A, 315 (1974). In the planar context, see C. 
Duval and P. A. Horvathy, Phys. Lett. B594, 402 (2004); Theor. Math. 
Phys. 144, 899 (2005). 

[12] J. Lukierski, P. C. Stichel, W. J. Zakrzewski, Annals of Physics (N. 
Y.) 260, 224 (1997); J. Lukierski, PC. Stichel and W.J. Zakrzewski, 
Ann. Phys. (NY.) 306, 78 (2003) 

[13] C. Duval and P. A. Horvathy, Phys. Lett. B 479, 284 (2000) ; C. Duval 
and P. A. Horvathy, Journ. Phys. A 34, 10097 (2001) 

[14] A. A. Deriglazov, JEEP 03, 021 (2003) 

[15] M. S. Plyushchay, Phys. Lett. B273, 250 (1991). M. S. Plyushchay, 
Phys. Lett. B 262 (1991) 71; P. A. Horvathy and M. S. Plyushchay : 
JHEP 06, 033 (2002); P. A. Horvathy, M. S. Plyushchay and M. Valen- 
zuela, Ann. Phys. 325 1931-1975 (2010) larXiv: 1001 .02741 

[16] J.-M. Souriau, Structure des systemes dynamiques, Dunod,Paris (1970); 

[17] E. Cartan, Legons sur les invariants integraux. Hermann: Paris (1922). 



14 



[18] L. D. Faddeev and R. Jackiw, Phys. Rev. Lett. 60, 1692 (1988). 

[19] P. A. Horvathy, Ann. Phys. 299, 128 (2002) 

[20] P. A. Horvathy, Journ. Math. Phys. 20, 49 (1979) 

[21] D. R. Grigore, Int. J. Mod. Phys. A 7, 7153 (1992) 

[22] G. Dunne R. Jackiw and C. A. Trugenberger, Phys. Rev. D41, 661 

(1990) ; 

[23] R. B. Laughlin,Phys. Rev. Lett. 50, 1395 (1983). 

[24] P. A. Horvathy, L. Martina, P. Stichel in " Noncommutative Spaces and 
Fields", P. Aschieri et al. ed.s , SIGMA 6 (2010), 060. 

[25] D. R. Grigore, Journ. Math. Phys.. 37, 460 (1996). 

[26] A. Ballesteros, M. Gadella, and M. del 01mo,Journ. Math. Phys. 
33, 3379 (1992); Y. Brihaye, C. Gonera, S. Giller and P. Kosihski, 
[hep-th/9503046 (unpublished)]. 

[27] C. Duval, Z. Horvath, P. A. Horvathy, L. Martina and P. Stichel, Mod. 
Phys. Lett. B20, 373 (2006) ; C. Duval, Z. Horvath, P. A. Horvathy, L. 
Martina and P. Stichel, Phys. Rev. Lett 96, 099701 (2006). 

[28] D. Xiao, J. R. Shi and Q. Niu, Phys. Rev. Lett. 95 (2005) 137204 

[29] A. Berard and H. Mohrbach, Phys. Rev. D 69, 127701 (2004); P. Gos- 
selin, F. Menas, A. Berard and H. Mohrbach, Europhys.Lett. 76 (2006) 
651-656. 

[30] P. A. Horvathy, Phys. Lett. A 359 (2006), 705-706 
[31] K. Yu. Bliokh, Phys. Lett. A351 123 (2006) ; 
[32] Snyder H.S., Phys. Rev. 71 (1947), 38. 

[33] R. Jackiw and V. P. Nair, Phys. Rev. D43, 1933 

(1991) M. S. Plyushchay, Phys. Lett. B248, 107 (1990); 

[34] R. Banerjee Eur. Phys. J. C 47, 541545 (2006); Phys. Rev. D 75, 
045008 (2006). 



15 



